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In the article of V.V. Anisovich et al., Yad. Fiz. 63, 1489 (2000), the Tf-matrix solutions for 
the wave IJ^'^ = 00++ were obtained in the mass region 450-1900 MeV where four resonances 
/o(980), /o(1300), /o(1500), /o(1750) and the broad state /o(1530 +250) are located. Based on these 
solutions, we determine partial widths for scalar-isoscalar states decaying into the channels tttt, KK, 
rjT), 1777', TTTTTTTr and corresponding decay couplings. 



- - - I. INTRODUCTION 

I In the combined if-matrix analysis was performed for the meson partial waves / J^*^ = 00++, 10++, 02++, 12++ 

. on the basis of GAMS data on Tr^p n^Tr^n, ijrin, rjri'n BNL data on -K^p — > KKn |^ and Crystal Barrel data 

' on pp{at rest) — > tt^tt^tt^, t^^Wj t^^t^^V 0- The positions of amplitude poles (physical resonances) were determined 

d . together with the positions of the /C-matrix poles (bare states) and bare-state couplings to the two-meson channels. 

^! The nonet classification of the bare qq states was suggested, and the possibilities for the location of the lightest scalar 

, glueball in the mass region 1200-1700 MeV were discussed. 

^^O > The if -matrix technique has an advantage of taking account of the unitarity condition constraints and, in this way, 

' of a correct incorporation of threshold singulatities into the scattering amplitude. Since the search for resonances is 

^ . always related to the investigation of analytical structure of the amplitude in the complex-mass plane by using data 

QQ ' at real masses, it is important to perform analytical continuation of the amplitude into the lower half-plane of the 

, complex mass, with correctly taken singularities on the real axis. 

However, the if-matrix amplitude does not include resonance parameters in an explicit form, so additional calcu- 

CN ] lations are needed to determine masses and couplings of real resonances. In paper [|l| the pole positions have been 

' found for the considered partial wave amplitudes (i.e. masses and total widths of resonances were found), but more 

\^ complicated calculation of couplings has not been done yet. The decay coupling constants are to be determined as 
residues of the pole singularities of the multi-channel amplitude. In the present paper we overcome this deficiency and 

' calculate coupling constants to the channels tttt, KK, r]r], 7777' and ttttttt: for the resonances /o(980), /o(1300), /o(1500), 

O-i, /o(1750), /o(1530 ^2^0)- This procedure provided us with partial widths for the decays of these resonances. The choice 

O •■ of the scalar-isoscalar sector for primary study follows from the interest in pursuing the destiny of the lightest scalar 

, glueball after its mixing with neighbouring states; one needs a knowledge of the couplings /o tttt, KK, ijr], -q-q' for 

■ all resonances over the mass region 1000-1800 MeV. 

^ \ The paper is organized as follows. 

In Section 2 the coupling constants are presented for the decays /o — > tttt, if if, 7777, 7777', TTTTTTTr, and partial widths 
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for the mesons /o(1300), /o(1500), /o(1750) and /o(1530i^°o) are determined. 

In Section 3 the resonance /o(980) is considered in detail: the results of the analysis |^ tell us that standard 
formulae for the description of resonances, such as Breit-Wigner or Flatte ones, in the case of /o(980) are unable to 
give simultaneously the values of the decay coupling constants and the position of the amplitude pole. We suggest an 
alternative form of resonance amplitude for /o(980) in which the important role is played by the prompt transition 
TTTT KK. In this Section another low-mass state, namely, the cr-meson, is also discussed. 
The results are summarized in the Conclusion. 



II. DECAY COUPLINGS AND PARTIAL DECAY WIDTHS 



We determine the coupling constants and partial decay widths using the following procedure. The 00++-amplitude 
for the transition a — s- 6, 

Aa~,b{s) , a, b = ttt:, KK,rir],r]ri' , TTTTTTTr (1) 

is considered as a function of the invariant energy squared s in the complex-s plane near the pole related to the 
resonance n. In the vicinity of the pole the amplitude reads: 

(«) («) 

A,^,{.s) = ^-\^e<' (2) 
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Here /i„ is the resonance complex mass /i„ = M„ — ir„/2 ; ga"-* and g^"-* are the couphngs for the transitions /o — > a 

and /o b. The factor exp (iO^ab) is due to a backgromid contribution which can be the non-resonance terms or tails 
of neighbouring resonances. We also write down in (H) the non-pole background term Bab- 

The partial width for the decay /o — ^ a is determined as a product of the coupling constant squared, gi"'' ^, and 
phase space, Pa{s), averaged over resonance density: 

r ( \ [ — gi"^ ^ Pa{s) 

s>sth 

Following 0] , we write down the phase space factor as follows: 

= ^ , (4) 

where ka is the relative momentum of mesons in the decay channel (for example, for the tttt channel p^^[s) = 
yj{s — Am'^)/s ). For the innnr channel the phase space factor was chosen in |l] to be the same as for the two- 
p-meson state at s < 1 GeV^ or be equal to 1 at s > 1 GeV^. The integration over s in (||) is carried out in 
the region above the a-channel threshold, s > Sth (for the tttt channel it is s > 4rn^). The resonance density 
factor, ({Re fx^ — s)^ + (Im/i^J^) guarantees rapid convergence of the integral (||). The normalization constant 
Cn is determined by the requirement that the sum of all hadronic partial widths is equal to the total width of the 
resonance: 

r(n) = ^r,(n). (5) 

a 

In paper three solutions for the wave IJ^'~" = 00++ have been found; they are labelled as I, II- 1 and II-2 (see 
Tables 4 and 5 in In practice Solutions II-l and II-2 give the same physical parameters of resonances, though 
they differ from parameters found for the if-matrix elements. In particular, in Solution II-2 the state /o°'''^(1600) may 
be identified as a gluonium, for the decay couplings satisfy all the requirements inherent to giueball state; in Solution 
II-l such a state is /o*'"^^(1230). For Solution I, the same bare state, /q°'''^(1230), should be considered as a gluonium. 

In Table 1 we show the values of partial widths for the resonances /o(980), /o(1300), /o(1500), /o(1750) and broad 
state /o(1530 ^250)- Partial widths for /o(1300), /o(1500), /o(1750), /o(1530 ^2^0) are calculated within standard 
formulae for the Breit-Wigner resonances (^), (^ and (||). The resonance /o(980) being located near the strong KK 
threshold needs a special consideration that is presented below. 

The decay coupling constants squared, gi""* ^, are shown in Table 2 for a = tttt, KK, rjrj, -qri' , Trmrn. The couplings 
are determined with the normalization of the amplitude used in ^ : for example, we write the tttt scattering amplitude 
(^ as ATTTr^Trir (s) = (^70 ^xp (2i(5o) — l)/2ip^7r(s), where 77° and 6q are the inelasticity parameter and phase shift for 
the 00++ TTTT-wave. The coupling constants gi"'' are found by calculating the residues of the amplitudes tttt — > tttt, 
KK, rjrj, rjr]' , ttttttt:. Also we check the factorization property for the pole terms by calculating residues for other 
reactions, such as KK KK. 

The position of poles in the complex- Af plane (M = y/s) is illustrated by Fig. 1. The complex-M area, where 
the JC-matrix fit Q may reliably reproduce analyticity of the amplitude, is inside a semi-circle depicted by dashed 
line. The poles which are a subject of the iC-matrix analysis and correspond to /o(980), /o(1300), /o(1500), /o(1750), 
/o(1530 ^2^0) located on the 3rd, 4th, Sth and 6th sheets of the complex-M plane. The resonance /o(980) is 
located near the strong KK threshold, therefore two poles are related to /o(980): the nearest one is on the 3rd sheet 
(M ~ 1014 - i39 MeV) and a remote pole on the 4th sheet (M ~ 936 - i238 MeV). Couphng constants for /o(980) are 
determined as residues of the nearest pole (on the 3rd sheet). The ijr]' threshold is weak for /o(1500), and because of 
that positions of poles on the Sth and 6th sheets are practically the same (note that couplings related to these poles 
nearly coincide). 

The X-matrix fit [Q| has been carried out in a broad mass interval, 450 MeV < M < 1900 MeV. This very fact 
allows us to believe that we deal with successfully reconstructed analytical amplitude which goes rather deeply into 
lower half-plane M, and this area is restricted by dashed line in Fig. 1. It is of crucial importance that the pole of 
the broad state /o(1530+^n) is inside this area, for this broad state plays a key role in a mixing of gg-mesons with 
the lightest giueball, see |^,^ for details. 
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III. THE LOW-MASS MESONS: Fo(980) AND a 
The two low-mass mesons, /o(980) and cr- meson, need special consideration and comments. 

The analysis [Q shows us that /o(980) cannot be described either by standard Breit-Wigner formula or its modifi- 
cation for the case of the nearly located strong KK threshold, that is, Flatte's formula j^]. Here we suggest another 
resonance formula for /o(980) which agrees with the results of 

In the compilation of Particle Data Group 1^, the cr-meson is denoted as /o(400 — 1200) that reflects a cumulative 
result obtained in a number of papers where the mass of a was found in this region or even higher. However, the 
analysis |0] definitely demonstrates the absence of poles in the 00++-amplitude at (600 < Re M < 1200) MeV, with 
an exception of poles for /o(980) - we will discuss the situation with cr-meson in this Section later on. 



A. Description of /o(980) 

For /o(980), the X-matrix fit |] gives us the position of the pole and coupling constant values, see Tables 1 and 
2. These parameters are sufficient to reconstruct the Breit-Wigner resonance amplitude. However, in case of /o(980) 
there exists a strong KK threshold near the pole, so the resonance term in the amplitude (^) should be suggested not 
as the Breit-Wigner pole but in a more complicated form. For the tttt — > tttt and KK — > KK transitions near /o(980), 

the following resonance amplitudes can be written instead of the Breit-Wigner pole term i?!"^' — 5a"''.9["V(/^n ~ 
entering equation (^): 



where 

F^2GGj,jif + f{ml~s) , 

D = ml^.s-tG' -i^^^^(Gl^ + iF) . (7) 

mo \ 

Here toq is the input mass of /o(980), G and are coupling constants to pion channels (tttt -|- tttttttt) and KK. 

The dimensionless constant / stands for the prompt transition KK tttt: the value f /mo is the "transition length" 
which is analogous to the scattering length of the low-energy hadronic interaction. The constants niQ, G, Gj^f^, f are 
parameters which are to be chosen to reproduce the /o(980) charactristics (position of pole s ~ (1.02 — i0.08) GeV^ 
and couplings to the channels tttt and KK, g^^ ~ 0.076 GeV^ and ~ 0.184 GeV^, see Tables 1 and 2). 
The TTTT scattering amplitude in the /o(980) region is now defined as 

A - p*^ nlTrTr.^^) , if • ^ (o\ 

^TTTV^TTTV — & ^ fo(9m) + 6 2 Sm - . [V,) 

This formula may be compared with equation (^): the background term in (H) is fixed by the requirement that tttt 
scattering amplitude below KK threshold has the form exp (i5) sin 5. 

At / the resonance amplitudes (||) turn into Flatte's formula 0, which is used rather often for the description 
of /o(980) . Still, it happened that position of the pole (complex mass value) as well as the amplitude residue in the 
pole, which have been determined in |^ and shown in Tables 1, 2, do not obey the Flatte formula but require / 7^ 0. 

We obtained two sets of parameters, with sufficiently correct values of the /o(980) pole position and couplings. 
They are equal (in GeV units) to: 

Solution A : rriQ = 1.000 , / = 0.516 , G = 0.386 , G^^ = 0.447 , 

Solution B: mo = 0.952 ,/ = -0.478 , G = 0.257 , G^^ = 0.388 . (9) 

The above parameters provide us with a reasonable description of the tttt scattering amplitude. The phase shift Sq and 
inelastisity parameter rjQ are shown in Fig. 2; the angle 9 for the background term in Solutions A and B, determined 
as 

= 01 + (— - 1)02 , (10) 
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is numerically equal to 



Solution A: Oi = 189°, 62 = 146° , 

Solutions: 6*1 = 147°, 02 = 57°. (11) 

Solutions A and B give significantly different predictions for 77°; however, the existing data j|,|l^ do not allow us to 
discriminate between them. 

Partial widths of /o(980) are calculated with the expression similar to (3) — with the replacement of the integrand 
denominator as follows: 

(Re/i2-s)2 + (ImMD'^l^l'- (12) 

For both sets of parameters (^) the calculated partial widths are close to each other. For example, using Solution II-2 
and the ^-set of parameters we have F^^r = 62 MeV, Fj^^ = 14 MeV, while for solution n-2 and the i?-set of the 
parameters one has Ftttt = 66 MeV, F^^ = 10 MeV. Fhe values of partial widths for /o(980) averaged over Solutions 
A and B are presented in Fable 1. 

Fhe total hadron width of /o(980) is defined in the same way as for the other /g-mesons, namely, by using the 
position of pole in the complex-M plane: the imaginary part of the mass is equal to a half-width of the resonance. For 
the Breit-Wigner resonance this definition is in accordance with what is observed from resonance spectrum (provided 
there is no interference with the background). If the resonance is located in the vicinity of a strong threshold, the 
observed resonance width can differ significantly from what is given by the pole position. In Fig. 3 one can see the 

magnitudes l^j^^^^ggj)' Pp7r(s), l-R/^(g^'o)^^' PP-R'('S) and gl.^rtot/\D{s)\^ for the parameter sets A and B: the width of 
peaks does vary, being less than the value determined by the complex mass of the resonance. 

Fhe only objective characteristic of the total hadron width is the position of pole in the complex-M plane, due to 
this reason we employ such a definition of total hadron width. Multiple variations of total width in the compilation 
1^ are just due to the absence of a proper definition of Ttot for resonances near the strong threshold. 



B. The light ct- meson 



Fhe light cr-meson reveals itself as a pole on the 2nd sheet: it is shown in Fig. 1 at M = (431 — i325) MeV (or, 
in terms of s which is more appropriate variable for light particles, at s = (4 — ilA)m'^ ) that corresponds to the 
magnitude obtained in Although this pole does not appear in the area of complex M, where the if-matrix fit 
101 reconstructs the amplitude rather reliably, it still deserves detailed comments. 

Fhe situation with cr-meson is as follows. The if-matrix representation allows us to reconstruct correctly the 
analytical structure of the partial amplitude in the physical region, at s > 4mJ by taking account of the threshold 
and pole singularities. The singularities related to forces (or left singularities, at s < 0), are not included directly 
into the iC-matrix machinery. This does not allow us to be quite sure about the results of the if-matrix approach at 
s < 4m^. Concerning the low-mass region, s > 4:m'^, an important result of the if-matrix fit is the absence of the 
pole singularity in the 00+"'" amplitude at 500-800 MeV. Here the TTTr-scattering phase i5q increases smoothly reaching 
90° at 800-900 MeV. A straightforward explanation of such a behaviour of i5q could consist in the existence of a 
broad resonance, with the mass about 600-900 MeV and width ^800 MeV (for example, see discussion in p2|jl^ ] and 
references therein). However, as was stressed above, the X-matrix amplitude does not contain pole singularities at 
500 < Re M < 900 MeV: the iiT-matrix amplitude has a low-mass pole only, which is located near the tttt threshold or 
below it. In |0] the presence of the pole near the tttt threshold was not emphasized, since the if-matrix solution does 
not guarantee a reliable reconstruction of the amplitude at s ~ 4m^. In jTll , in order to restore analytical structure 
at s '--^ 4m^, the left-hand-side singularities were accounted for on the basis of the dispersion relation A^/Z?-method. 
The TTTT scattering N/ Z?-amplitude was represented at M < 900 MeV being sewed with the i^-matrix solution Q at 
450 < M < 900 MeV. The 7V/_D-amplitude reconstructed in this way has a pole near the tttt threshold, thus proving 
that qualitatively the results of ^ are also valid for the region s ~ 4m J. The pole of the A^/_D-amplitude is 
shown in Fig. 1. 

It is worth mentioning that the low-mass location of the tr-meson pole was also obtained in a set of papers, where the 
low-energy tttt amplitude has been investigated by taking into account the left-hand cut as a set of meson exchanges. 
These papers include: (i) dispersion relation approach, s ~ (0.2 — i22.5)m^ p^ , (ii) meson exchange models, s ~ 
(pi) - il7.8)ml [[|, s =i (p.5 - il3.2)ml [|l6|, s ~ (2.9 - ill.8)ml 0, (iii) linear cr-model, s ~ (2.0 - il5.5)ml 
At the same time, in p9|-^ the pole position was found in the region of higher s, at s > 7m^, that reflects the 
ambiguities of approaches which treat the left-hand cut as a known quantity. 
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As to finding out the location of cr-meson on the basis of the available experimental data, one should make general 
remark. Since the width of the cr-meson is rather large, it is necessary to fit to data in the energy interval which 
is much larger than a total width of the a-meson. For example, to speak about cr-meson with a mass ~ 900 
MeV and half- width r/2 ~ 400 MeV, one should fit to data in the interval 300 MeV < M < 1400 MeV and at the 
same time to take a correct account of the nearest singularities, which are poles corresponding to /o(980), /o(1300) 
and presumably /o(1500) as well as threshold singularities tttt, KK, rj^rj' and mrmr. Concerning nnmr, one should 
have in mind that the contribution of this channel is significant starting from 1300 MeV, so this channel is absolutely 
necessary. Such demands towards the fit of experimental data have been fulfilled in no paper under discussion, with 
an exception for p]JlO[|. 

IV. CONCLUSION 

We have obtained partial decay widths for five scalar-isoscalar states /o(980), /o(1300), /o(1500), /o(1750), 
/o( 1530^250) t>y calculating the decay couplings as residues of pole singularities in the if-matrix amplitude [Q: 
positions of poles in the complex-M plane are shown in Fig. 1. The pole which corresponds to the light ct- meson is 
also shown in Fig. 1: it was not included into the JsT-matrix calculation procedure directly, being close to the left-hand 
cut; the discussion of its status can be found in and references therein. 

The results of our calculations of partial decay widths are presented below (the magnitudes are given in MeV units): 



Resonance 


r 




i T]7J 




r 


rtot/2 


/o(980) : 


64 ±8 


12 ± 1 








4±2 


40 ±5 


/o(1300) : 


46± 12 


5±3 


4±2 






171 ± 10 


113 ± 10 


/o(1500) : 


37 ±2 


7±2 


4± 1 


0.2 ± 


0.1 


79 ±6 


62 ±3 


/o(1750) : 


74+15 


11+17 


7± 1 


3± 


1 




qo+20 


/o(1530tf5o) 


380 ± 15 


185 ± 10 


40 ±5 


1 ± 


1 


560t?^° 


ccnn+120 

O9U_200 



(13) 

The values shown for partial widths as well as decay coupling constants of Table 2 need some comments. 

The comparison of the hadron decays /o(980) — > KK and /o(980) — > tttt points to a large ss component in /o(980). 
The analysis of radiative decays 0(1020) — > 7/o(980) and /o(980) — > 77 shows also that the ss component in 
/o(980) is large: with the /o(980) flavour wave function written as nficos ip + ss sin ip, the radiative decay widths give 
either (p ~ —48° or ~ 86° (solution with negative (p is more preferable). When the decay processes are switched off, 
/o(980) transforms into /q°'''^(720± 100), with ipbare — —70° (corresponding pole trajectory in the complex-Af plane 
is shown in Fig. 10 of 0|). We see that the decay processes and related change of the state do not diminish the ss 
component strongly. 

An opposite situation takes place with /o(1750). After switching off the decay channels, this resonance transforms 
into /q'"''^(1810 ± 30) which is dominantly ss: (pbare — 90° for Solution I and (fitare — —60° for Solution II. However, 
partial decay widths of /o (1750) (or decay coupling constants given in Table 2) unambiguously prove that ss component 
in /o(1750) decreased strongly due to a mixing with other states after the onset of the decay processes. It is possible 
to guess that this ss component has flown into the broad state /o(1530t25o)- tti^ ratio r^^/r^^r for /o(1530t25o) 
does not contradict such an assumption. Such a scenario looks rather intriguing, in particular when taking account 
of the fact that the broad state /o(1530jl25o)' according to [|l|, is a descendant of a pure glueball (see also p, ^ , ^ ) . 
However, the study of the mixing of qq-state with the glueball is beyond the frame of this article; it will be investigated 
elsewhere. 

For /o(980), the obtained magnitudes for the complex mass and decay couplings and gj^, demonstrate a failure 
of the Flatte formula. We suggest an alternative description of /o(980) which explores, as an addition to the pole 
term, the amplitude for the prompt transition tttt — ^ KK. 
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TABLE I. Partial widths of scalar-isoscalar resonances (in MeV units) in hadronic channels tttt, KK, rjri, rirj' and mrmr for 
different -matrix solutions of Ref. [1]. 





TTTT 


KK 


VV 


vv' 


TTTT TTTT 


pole position 


solution 


/o(980) 


71 


13 


- 


- 


6 


1006 - 


i 45 


I 




56 


10 


- 


- 


2 


1020 - 


i 34 


II- 1 




64 


12 


- 


- 


3 


1015 - 


i 39.5 


II-2 


/o(1300) 


71 


12 


7 


- 


170 


1307- 


i 130 


I 




34 


2 


2 


- 


171 


1296 - 


i 104 


II- 1 




34 


2 


2 


- 


171 


1296 - 


i 104 


II-2 


/o(1500) 


38 


9 


4 


0.3 


83 


1495 - 


i 67 


I 




37 


6 


4 


0.1 


83 


1498 - 


i 60 


II- 1 




37 


6 


4 


0.1 


72 


1495 - 


i 60 


II-2 


/o(1750) 


45 


28 


6 


4 


29 


1781 - 


i 56 


I 




88 


2 


7 


3 


120 


1817 - 


i 110 


II- 1 




90 


2 


7 


3 


123 


1817 - 


i 113 


II-2 




376 


173 


42 


2 


818 


1609 - 


i 706 


I 




377 


196 


39 


1 


434 


1427 - 


i 526 


II- 1 




376 


196 


39 


1 


438 


1430 - 


i 525 


II-2 



TABLE 11. Coupling constants squared (in GeV^ units) of scalar-isoscalar resonances to hadronic channels tttt, KK, r)r}, r)r}' 
and nnmr for different iiT-matrix solutions of Ref. [1]. 





TTTT 


KK 


VV 


vv' 


TTTT TTTT 


solution 


/o(980) 


0.076 


0.180 


0.075 




0.009 


I 




0.076 


0.186 


0.072 




0.004 


II- 1 




0.076 


0.186 


0.072 




0.004 


II-2 


/o(1300) 


0.050 


0.015 


0.012 




0.124 


I 




0.026 


0.002 


0.003 




0.132 


II- 1 




l).()2() 


0.002 


0.003 




().l:!2 


11-2 


/o(1500) 


0.032 


0.010 


0.005 


0.012 


0.070 


I 




0.038 


0.009 


0.007 


0.006 


0.074 


II- 1 




0.038 


0.009 


0.007 


0.006 


0.074 


II-2 


/o(1750) 


0.039 


0.029 


0.007 


0.030 


0.025 


I 




0.086 


0.003 


0.009 


0.028 


0.117 


II- 1 




0.086 


0.003 


0.009 


0.028 


0.117 


II-2 


/o(1530+^^o) 


0.329 


0.229 


0.061 


0.022 


0.764 


I 




0.304 


0.271 


0.062 


0.016 


0.382 


II- 1 




0.304 


0.271 


0.062 


0.016 


0.382 


II-2 
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FIG. 1. Pole positions (full circles) in the complex- Af plane {M = ^/s). Solid lines stand for cuts related to the threshold 
singularities {nn, nnm:, KK, rjrj and rjrj'). Two poles, which correspond to /o(980), are shown: on the 3rd and 4th sheets. 
On the 5th sheet the poles for /o(1300), /o(1500) and broad state /o(1530l25o) ^'^'^ located (for the broad state the pole 
stands at (1430 — i525) MeV, that is the mass for Solution II of the fit [1]). On the 6th sheet there is a pole for /o(1750). 
The dashed semi-circle restricts the area where the matrix fit [1], which was carried out on the real axis in the interval 
450 MeV < M < 1900 MeV, can give a reliable reconstruction of analytical amplitude. 
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FIG. 2. Reaction tttt — » tttt: description of 5q and rjo in the region of /o(980). Solid and dashed curves correspond to the 
parameter sets A and B. Data are taken from [9] (full squares) and [10] (open circles). 
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FIG. 3. a,b) The magnitudes |.R/^(^8o)VP'r('5) (solid curve) and Pk{s) (dashed curve) for the parameter sets 

A and B. Visible peaks in the tttt spectra have the total widths ~ 60 MeV (set A) and ~ 45 MeV (set B). c,d) Values 
gi.^Vtotl\D(ay\' (where g^^Vtot ~ 0.006 GeV^) for the parameter sets A and B. Visible total widths of the peaks are ~ 55 MeV 
(set A) and ~ 45 MeV (set B). 
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